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ABSTRACT 

Aims. The hadronic vs. leptonic origin of the gamma-ray emission from the Supernova Remnant RX J1713. 7-3946 is discussed 
both in the light of new observations and from a theoretical point of view. 

Methods. The existing good spatial correlation of the gamma-ray emission and the nonthermal X-ray emission is analyzed 
theoretically. In addition, the recently published new H.E.S.S. observations define the energy spectrum more precisely, in 
particular at the high and low energy ends of the instrument's dynamical range. There now exist much more constraining X-ray 
observations from Suzaku that extend substantially beyond 10 keV. These new data are compared with the authors' previous 
theoretical predictions, both for dominant hadronic and for simple inverse Compton models. 

Results. Apart from the well-known MHD correlation between magnetic field strength and plasma density variations, emphasized 
by the wind- bubble-structure of the remnant, it is argued that the regions of magnetic field amplification also are correlated 
with enhanced densities of accelerated nuclear particles and the associated streaming instabilities. Therefore a correlation of 
nonthermal X-ray and 7-ray emission is not only possible but even to be expected for a hadronic emission scenario. A leptonic 
origin of the gamma-ray emission would require an implausibly uniform strength of the magnetic field. The observational and 
theoretical inferences about substantial field amplification in this remnant agree very well with the recent X-ray and 7-ray 
observations. 

Conclusions. All this argues strongly for the dominance of hadronic 7-rays in the 7-ray emission spectrum and a fortiori for an 
overwhelming contribution of nuclear cosmic ray particles to the nonthermal energy in this remnant. 

Key words. (ISM:) cosmic rays - acceleration of particles - shock waves - supernovae individual(SNR RX J1713. 7-3946) - 
radiation mechanisms:non-thermal - gamma-rays :theory 



1. Introduction 

Hard X-ray emission from young Supernova Remnants 
(SNRs) was first detected in SN 1006 and interpreted 
as synchrotron emissi on from electrons in the energy 
range of tens of TeV ( Kovama et al. . 1995). Such emis- 
sion meanwhile has been found in many young shell- 
type objects. These electrons are also expected to pro- 
duce very high energy (VHE) 7-rays in the > 100 GeV 
range by inverse Compton (IC) collisions with low en- 
ergy photons from the ambient radiation field, in particu- 
lar the Cosmic Microwave Background (CMB). In regions 
of high gas density this may also be accompanied by a 
non-negligible emission of nonthermal Bremsstrahlung 7- 
rays It is therefore not a priori clear whether the VHE 
7-rays detected in several shell-type SNRs up to now 
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Enomoto et ail . 120061) . and RCW 86 (jHoppe et all . l2007f) 
- are indeed the result of inelastic collisions of energetic 
nuclei with gas atoms, as expected theoretically if a size- 
able fraction of the SN explosion energy goes into nuclear 
cosm ic rays (CRs) (jPrurv et al .], Il994t iBerezhko fc Volkl . 
19971 ). or whether they are merely due to IC collisions of 
the X-ray synchrotron emitting electrons. 

The latter inference seems at least in a rough sense pos- 
sible if the effective magnetic field strength in these objects 
is sufficiently low. Then the energy spectrum of acceler- 
ated electrons, calculated from the observed synchrotron 
spectrum, contains sufficiently many particles so that a 
strong enough IC 7-ray emission results. Since magnetic 
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fields are difficult to determine, a possible phenomeno- 
logical approach is to assume that the interior (post- 
shock) magnetic field in the SNR is the result of MHD- 
compression of a circumstellar field whose strength is at 
best a few /j,G. For RX J1713. 7-3946 and RX J0852.0-4622 
such a procedure works, at least as far as the overall ampli- 
tude of the 7-ray spectrum is concerned. For Cas A this 
is true as well. However, all these remnants presumably 
correspond to core collapse SN explosions into the wind 
bubble of a massive progenitor star, with a rather complex 
circumstellar magnetic field structure. In contrast, for the 
remnants of several type la SNe which could not be de- 
tected so far in 7-rays such as Tycho's SNR, SN 1006, and 
Kepler's SN, this scheme leads to large overpredictions 
for the 7-ray flux, unless the mean interior magnetic field 
strength B is assumed to be significantly higher than the 
value ~ 10 fiG expected from a ^ 5 /iG circumstellar field 
for a strong quasi-parallel, adiabatic shock with a com- 
pression ratio of 4 (Volk et al. 2007, 2008). Since electron 
acceleration cannot be responsible for such an amplifica- 
tion of the interstellar field, this makes a pure electron 
acceleration scenario generally untenable in our view. 

The acceleration of nuclear particles on the other 
hand cannot be deduced in a similar way from some 
other electromagnetic emission - even though the mag- 
netic field strength and the ion injection rate are in- 
ferred from the form o f the electron synchrotron spectrum 
(Berezhko et al. , 2002). The energetic nuclear component 
needs to be calculated theoretically, and the amplitude of 
the accelerated particle distribution can only result from 
a nonlinear theory. This has been done successfully in a 
number of cases and shows that the overall energy in ac- 
celerated nuclei is indeed a large fraction - about 10% - of 
the available hydrodynamic explosion energy, and that the 
local efficiency at those parts of the shock surface, where 
ion injection is effec tive, amounts to even about 50 % (see 
e.g. iBerezhkol 120051 for a review). 

Given this high acceleration efficiency, at rea- 
sonable thermal gas densities the expected rate of 
hadronic collisions with gas atoms and the resulting 7- 
ray production by 7r°-decay turns out to be consis- 
tent with the observations in the 7-ray dete c ted ob - 
jects analyzed so far: Cas A dBerezhko et all . 12003b ). 
RX J1713. 7-3946 dBerezhko fe VolkL l2006l| [hereafter re- 
ferred to as BV06]. dBerezhko fc Vdlkl . 120071). RX J0852.0- 



20031: iBerezhko et all l2003al iBerezhko fc Vol! l2004al : 



Volk et all. 



2005UParizot et alll2006l ). The field strengths 



4622 dBerezhko et all . 120081 )" Also for the other objects 
that we have analyzed the oretically - Tycho's SNR 



dVolk et al.l.l2002ll2005l.l2008h. S N 1006 dBerezhko et al 



200 21 iKsenofontov et all 120051). and Kepler's SNR 



dBerezhko et all 20061 ) - and for which only upper lim- 
its exist so far, the hadronic 7-ray emission is expected to 
dominate the IC fraction, even though for the low-density 
object SN 1006 only by a modest margin. 

Within the errors, the theoretically derived mag- 
netic field strengths also agree with those deduced from 
the observations of filamentary X-ray synchrotron fea- 
tures, often fou n d at t he outer rims of t he SNRs (e.g 
Vink fc Lanhngl . l2003t lLong et all 120031 : iBamba et al 



typically are an order of magnitude greater than the com- 
pressed circumstellar field, substantially reducing the IC 
spectrum amplitude for a given synchrotron spectrum and 
thus tending to yield a negligible leptonic contribution to 
the observed VHE 7-ray flux. 

Although the dominant acceleration of nuclear parti- 
cles is clearly favored by what has been said up to now, 
leptonic scenarios have been studied in so me detail (e.g. 
IPorter et alll2006HKatz fc Waxmanll2008l) . Especially for 
the experimentally best-studied SNR RX J1713. 7-3946 
the question is whether an IC spectrum, scaled in ampli- 
tude from the synchrotron spectrum in the sense described 
above, is compatible with the form of the observed 7-ray 
spectrum. A second question regards the observed spatial 
correlations of the morphology in hard X-rays, assumed to 
be the result of synchrotron radiation, and in 7-rays. Such 
a correlation has been observed in RX J1713. 7-3946 and 
might at first sight be attributed to a common leptonic 
population of energetic particles. It remains to be seen 
whether such an inference holds upon deeper scrutiny. 

Another, not quite resolved question concerns the 
absence of thermal X-ray emission in the two SNRs 
RX J1713. 7-3946 and RX J0852. 0-4622 that are spatially 
resolved at TeV energies. Also the radio synchrotron emis- 
sion is very weak in comparison to the strong X-ray syn- 
chrotron emission in b oth sources. Since the earliest obser- 
vations and analyses ( Slane et al. . 1999), the explanation 
for this situation has been the assumption that the SN ex- 
plosion occurred into the very rarefied stellar wind bubble 
of a massive progenitor star. The model of BV06 specifies 
the bubble structure in detail. It demonstrates that the 
main gas heating and particle acceleration occurs beyond 
the wind region and bubble, when the SNR shock propa- 
gates into the swept-up, radiatively cooled shell of inter- 
stellar matter that was generated by the forward shock 
driven into the ambient ISM by the expansion of the wind 
bubble. Therefore the wind bubble has two main effects: 
the SNR can reach a large size quickly and it is then rather 
quickly decelerated in the dense shell. During propagation 
of the SNR shock wave through the wind and bubble the 
gas density is very low, resulting in a very low thermal X- 
ray emission of this material. The shock's late encounter 
of the massive shell heats the nuclear particles, but leaves 
little time for the subsequent relaxation in the postshock 
region of the SNR that heats the thermal electrons by 
Coulomb collisions with the heavy ions. In this way a 
minimum of thermal X-ray emission from the remnant 
is combined with a maximum of hard X-ray synchrotron 
and hadronic 7-ray emission from the SNR. In addition, 
the shock is nonlincarly modified. For a SNR propagat- 
ing into a uniform medium with a uniform magnetic field 
this implies that the larger part of the shock surface cor- 
responds to a quasi-perpendicular sho ck with a s t rongl y 



reduced injection of nuclear particles ( Volk et al. . 2003f ). 



Suprathermal injection of ions is only possible in the quasi- 
parallel shock regions. If the spatial scales of the quasi- 
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perpendicular regions are large enough, the cross-field dif- 
fusion of the highest-energy particles from the neighbor- 
ing quasi-parallel shock regions also does not reach deeply 
into these quasi-perpendicular regions. In the correspond- 
ing magnetic flux tubes no magnetic field amplification oc- 
curs either and the shock remains essentially un-modified 
there. This means that in the quasi-perpendicular regions 
the shock dissipation and therefore the gas heating oc- 
curs in a locally unmodified shock with the overall shock 
speed, leading to a rather high gas temperature. In the 
case where a radiatively cooling shell of a wind bubble is 
the main obstacle for the SNR expansion, the situation 
may be different. The MHD instabilities of such a shell 
probably break it into many small regions with strongly 
varying field directions. Then the spatial scales separating 
the quasi-perpendicular from the quasi-parallel shock re- 
gions may become small enough that crossfield diffusion 
can smear out the quasi-perpendicular regions and par- 
ticle accelerati on occurs practic a lly everywh ere over the 



shock surface (jVolk et all . l2007t IVolkl . 120081 ). In the ex 



treme this implies shock modification over the entire shock 
region and thus a reduced gas heating due to the subshock 
dissipation onljQ. Compared to an un-modified shock of 
the same overall speed, the gas temperature is then dimin- 
ished by a factor « of /a -2 ~ 0.25, where a and a s denote 
the overall shock compression ratio and the subshock com- 
pression ratio, respectively. This reduces the emission of 
soft thermal X-rays drastically. 

Such a configuration is not easily analysed in detail 
with standard methods of X-ray astronomy. However, this 
also leaves uncertainties in the evaluation of the model's 
validity. Here, as in BV06, we proceed under the assump- 
tion that the concrete wind bubble model, or its eventual 
improvement, is consistent with the fact that no thermal 
emission has been found up to now. 

The purpose of the present paper is a discussion of the 
above two questions for RX J1713. 7-3946 , even though 
we expect analogous arguments to hold for the other ob- 
jects mentioned. In section 2 we discuss the correlations 
between the gas density and the magnetic field strength 
that are possible, or are even to be expected, in the cir- 
cumstellar medium before the SN explosion. In an approx- 
imate way we also show how variations in gas density lead 
to correlated variations in nuclear CR pressure and in the 
hadronic 7-ray emission. However, we shall argue that the 
magnetic field direction influences particle injection into 
the shock acceleration process, so that circumstellar den- 
sity enhancements are only a necessary but not a sufficient 
condition for enhancements of the 7-ray emission. Under 
most circumstances electron synchrotron emission is then 
expected to increase with hadronic 7-ray emission as well. 
We also briefly discuss the alternative purely leptonic sce- 
nario. In section 3 we compare our calculations for the 
hadronic 7-ray emission, the synchrotron emission, and 
for the IC and Bremsstrahlung 7-ray emission with the 



latest X-ray and 7-ray measurements for RX JI7I3. 7-3946 
. The results are finally compared with a purely leptonic 
scenario, calculated by assuming ion injection and thus 
also magnetic field amplification to be negligibly small. 
It turns out that a hadronic origin of the 7-ray emission 
is consistent with all measurements and with our theory, 
whereas the leptonic scenario runs into serious difficulties. 
Section 4 contains our conclusions. 

Very recently an in d epend ent discussion of these ques- 
tions by iTanaka et all (|2008l) has come to our attention 
Their paper is more observation oriented and is partly 
complementary to ours. It uses phenomenological esti- 
mates for the accelerated particle spectra. However their 
conclusions are similar to ours. 

2. Expected morphological correlations 

In this section we investigate the spatial correlations that 
arise from compressions or expansions of the thermal gas 
with a frozen-in magnetic field, in a given radiation field. 
In addition we argue that variations of the nuclear ener- 
getic particle density - which are most likely the result of 
spatial variations of the injection rate into the accelera- 
tion process - lead to positively correlated enhancements 
of the magnetic field strength as a result of magnetic field 
amplification by CR streaming instabilities. 

2.1. Compressions/de-compressions of the thermal 
plasma 

Consider briefly quasi-static compressions / de- 
compressions of the gas density by plasma motions 
across the magnetic field in the MHD limit. 

In this case the conservation of magnetic flux com- 
presses the field together with the thermal gas, increasing 
at the same time the target for electrons to produce syn- 
chrotron emission as well as the target for energetic nuclei 
in inelastic collisions with the gas atoms. In both cases 
the densities of the X-ray synchrotron emitting electrons 
and of the 7-ray emitting nuclei vary approximately in 
the same way, because the electrons and nuclei concerned 
have roughly the same energy. Therefore their spatial dis- 
tributions are essentially the same. 

In the I-fluid MHD-approximation the mass density p 
and the magnetic field vector B are related to the mass 
velocity u through the conditions of conservation of mass 
and magnetic flux: 



dp 

at 



u ■ Vp 



dp 
dt 



-pV ■ u 



OB 

— = V x (u x B). 

dt x ' 



(1) 



(2) 



1 The given 7-ray flux requires in addition a reduction of the 
gas density which reduces the thermal emission. 



Being interested primarily in the variations of p and the 
magnitude B of B, it is more useful for our considerations 
to reduce Eqs.(l-2) to equations for B 2 and B/p alone, 
using V • B = 0: 

1 d R 2 

2«^- = -V.«-B.(B.V)« (3) 
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dln(B/p) 



= B~ 2 B ■ (B ■ V)it, 



(4) 



where the r.h.s. of Eq.(3) is nonzero only for an expan- 
sion/compression perpendicular to the direction of B and 
Eq.(4) shows that the ratio B/p changes only in an ex- 
pansion/compression parallel to B. 

With these relations in mind we consider the circum- 
stellar medium of a single massive star before it is reached 
by the SNR shock. The stellar wind produces an expand- 
ing hot, low density bubble. This dynamical evolution 
compresses the ambient interstellar medium (ISM) in an 
outer shock wave in an approximately spherically symmet- 
ric manner. It is clear that the compression of this external 
medium will create a shocked shell with the gas density 
increasing radially outwards and, over most of the solid 
angle (where the interstellar B- field is not radial), also an 
increasing magnetic field strength. All pre-existing den- 
sity inhomogcneitics in the external medium in the form 
of clouds will by themselves have generally higher B-fields, 
wherever the density is high. The subsequent compres- 
sion by the wind bubble will tend to enhance these cor- 
related density and B-field contrasts. Azimuthal density 
variations will probably form upon this compression as 
well, by azimuthal instabilities of the radiatively cooling 
swept-up shell, regardless of pre-existing clouds. All this 
implies spatially correlated variations of p and B in the 
medium upstream of the SNR shock. 

As the theoretical model of BV06 suggests, the rem- 
nant of the subsequent SN explosion should be old enough 
that the swept-up mass of interstellar matter dominates 
the mass enclosed in the wind bubble as well as the ejected 
mass from the explosion. Then we have the situation that 
most of the emitting medium is a shock- modified ISM that 
originally had a structure as described above. Its nonther- 
mal emission dominates the 7-ray emission as well as the 
synchrotron emission of the SNR, as we observe it today. 

2.2. Injection and acceleration of nuclear particles at 
the SNR shock 

If the magnetic field lines in front of the accelerating for- 
ward SNR blast wave are quasi-perpendicular to the shock 
normal, suprathermal ions cannot escape from the down- 
stream region (i.e. from the remnant interior) into the up- 
stream region to start being diffusively accelerated. On 
the other hand, nuclear particles can be readily injected in 
quasi-parallel shock regions where the field lines are more 
or less parallel to the shock normal. And the number of 
ions injected into the acceleration from each gas volume is 
to a first approximation proportional to the thermal gas 
density. In such field regions a high energy density of accel- 
erated nuclear particles will be built up by the acceleration 
process. This particle population will try to escape into the 
upstream medium and thereby excite streaming instabil- 
ities which lead to magnetic field amplification in the as- 
sociated magnetic flux tubes. The enhanced field amounts 
at the same time to an increased target strength for syn- 



chrotron emission by electrons, while the enhanced ener- 
getic nuclear particle intensity will increase the number of 
inelastic collisions per unit volume at given gas density. 
As already mentioned in the Introduction an inhibition 
of acceleration in quasi-perpendicular shock regions may 
not occur in the majority of quasi-perpendicular regions 
of the shock because of cross- field diffusion. However, in 
principle this effect is there and may occur especially at 
pre-existing interstellar clouds encountered by the SNR 
shock. 

The injection of energetic electrons into the accelera- 
tion process is not very well understood, and is possibly 
not correlated with the ion injection. Electron injection 
might actually occur everywhere over the SNR shock sur- 
face. It may or may not increase with the plasma density. 

Let us for a moment assume that electron injection 
increases strictly monotonically with the upstream gas 
density. Then we would expect a nonlinear correlation 
between the hard X-ray and the 7-ray emissions in a 
hadronic model (see also section 2.3), since electrons do 
not play a role in the acceleration energetics and the re- 
sulting shock modification, acting de facto as test par- 
ticles. Therefore we expect the amplitude of the acceler- 
ated electron energy distribution to increase with p, even 
though its form is entirely determined by the accelerating 
ion component and is in fact equal to the form of the ion 
energy distribution. The reason is that only the nuclear 
component - with its dominant energy density - deter- 
mines the shock structure, and therefore the form of the 
accelerated energy spectrum of all ultrarelativistic parti- 
cles whose energy is high compared to the proton rest 
energy, including the electrons (as long as their radiation 
loss can be neglected). In fact, the lack of saturation of 
the electron acceleration will allow the energy density of 
nonthermal energetic electrons to increase strictly mono- 
tonically with their injection rate, i.e. with gas density, 
in contrast to the saturation behavior of the nonthermal 
nuclear particles which limits their injection rate. 

Even if the electrons were injected - and thus also 
accelerated - everywhere with the same intensity, then 
the amplification of the magnetic field by the ions in 
the quasi-parallel shock regions will automatically lead to 
an enhancement of the synchrotron emission in these re- 
gions, suggesting a correlation of the synchrotron emis- 
sion with the hadronic 7-ray emission. This should be 
qualitatively similar to the case in the bipolar remnant of 
SN 1006, even though probably less pronounced in wind 
bubbles. Th e bipolarity of SN 1006 has been discussed 



bubbles, ine bipo larity 01 bl \ lUUo lias been discussed 
in detail by Volk et all (|2003f ). For an experimental dis- 
cussion of the SN 1006 magnetic field configuration, see 



Rothenflug et all flgOOJJ) 



Summarizing the arguments so far, we might expect a 
good and perhaps even nonlinear spatial correlation be- 
tween the synchrotron emission and the production of en- 
ergetic nuclear particles, as a result of pre-existing corre- 
lated density and magnetic field strength inhomogeneities, 
as well as from field amplification in quasi-parallel shock 
regions. However, this last condition on nuclear particle in- 
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jection implies that not all gas density enhancements need 
to show up as regions of enhanced nuclear energeti c parti- 
cle d ensity and amplified magnetic field strength (jPlagal . 
20081) . In other words, not all clouds in the environment 
of RX J1713. 7-3946 need to be regions of enhanced syn- 
chrotron and hadronic 7-ray emission. Some might just 
be shielded against the injection of ions by an unfavorable 
magnetic field direction relative to the shock normal. 



2.3. Hadronic j-ray emission in the Quasi-Sedov phase 

The previous arguments did not consider the explicit mod- 
ulation effect of upstream density inhomogeneities in the 
upstream medium on the strength of nuclear particle ac- 
celeration and on the resulting hadronic 7-ray emission. 
To address this question we again make use of the as- 
sumption that the evolutionary phase of RX J1713.7- 
3946 is dominated by the swept-up mass from the wind 
bubble-structured interstellar medium. Then the SNR in- 
terior should be subsonic with a roughly uniform total 
gas and particle pressure Ptot oc piV^, as known from 
the Sedov solution in the case without CR acceleration. 
Let us call this phase the Quasi-Sedov phase. If for exam- 
ple, in a situation that deviates from spherical symmetry, 
the shock reaches regions of different gas density over its 
surface, then the local shock velocity V s there will vary 
with the local upstream gas density p\ as V s oc Pi°' 5 . In 
the Quasi-Sedov phase the nuclear energetic particle pres- 
sure P c is an only slowly varying function of time, locally 
equal to s everal 10 percent of Ptot if ion injection is effi- 
cient (e.g. lKsenofontov et all [2005). And this nonthermal 
energy density is primarily concentrated in the highest- 
energy particles of the spectrum, i.e. in the VHE range. 
This implies an essentially uniform VHE hadronic 7-ray 
emissivity q 1 oc P c in these regions and therefore a local 
hadronic 7-ray production rate Q 7 = q-yPd oc p±, since 
the downstream gas density pd is a fixed fraction of the 
upstream density p\ for a strong shock. 

Therefore the hadronic 7-ray emission is approxi- 
mately proportional to the local gas density. Unless the 
electron acceleration is anticorrelated with the upstream 
gas density, the X-ray synchrotron emission will be corre- 
lated with the hadronic 7-ray emission, because B and p 
are almost always spatially correlated. 

2.4. Spatial correlations in the leptonic scenario 

In the highly likely case that the leptonic 7-ray emission 
is dominated by the IC emission, the 7-ray morphology 
is basically determined by the spatial distribution of the 
radiating electrons. On the other hand, the synchrotron 
emission is proportional to the product of the energetic 
electron density and the local magnetic field energy. In 
the hypothetical case that the hadronic 7-ray emission is 
negligible, the observations would of course require a good 
spatial correlation between X-ray and IC 7-ray emission. 
And, the better this correlation in such a leptonic scenario, 



the more the magnetic field strength must therefore be 
spatially uniform. 

This is a most unlikely situation for RX J1713. 7-3946 
given the fact that the gas density appears quite variable 
around the SNR shell and the remnant appears to be inter- 
acting with at least so me of these molecular clouds, as CO 



obser vations suggest (jFukui et all 120031 : iMoriguchi et al 



2005). This makes a purely leptonic scenario already em- 



pirically highly unlikely, independently of any theoretical 
arguments. 



3. Comparison of the latest 7-ray spectrum with 
the previous theoretical model of BV06 

The latest version of the H.E.S.S. 7-ray spectrum 
( Aharonian et al. , 2007a) presented in Fig.l has not only 
increased statistical accuracy, but also a flat (hard) spec- 
tral shape at the lowest energies of 250 GeV, and a smooth 
extension and fall-off towards the 100 TeV region. For the 
evaluation of this spectrum the reflectivity changes of the 
H.E.S.S. mirrors have been taken into account, leading to 
a roughly 15 percent increase in flux compared to the 2005 
spectrum. 

Note that the theoretical 7-ray spectra, calculated in 
BV06 and also presented in Fig. 2, correspond to the as- 
sumptions that RX J1713. 7-3946 was a core collapse su- 
pernova SN of type Il/Ib with a massive progenitor and 
explosion energy E sn = 1.8 x 10 51 erg, that it has an age 
of w 1600 yr and is located at a distance of « 1 kpc. 
Although this general scenario corresponds to the conclu- 
sions of most other authors, significantly lar ger distances 



also have been considered in the literature (Slaneetal 



1999). For details, see the discussion in BV06. 



Since the theoretical model cannot well determine the 
spectral amplitude, for the reasons given in BV06, the 
above-mentioned flux increase is not relevant in a com- 
parison of theoretical and observational spectra^. 

However, the forms of the spectra agree remark- 
ably well. We note that the inferred leptonic IC and 
Nonthermal Bremsstrahlung spectra are depressed by 
more than an order of magnitude relative to the observed 
spectrum. They cannot explain the observations, if the 
magnetic field is indeed amplified to the degree assumed 
in the theory and supported by the upper limit for the 
thickness of the synchrotron filaments which one can de- 
riv e (BV06) from the X MM-observations of this remnant 
bv iHiraga et al.1 ( 20051 ). The resulting lower limit for the 
magnetic field strength of 65 p,G has recently been sup- 
ported by Chandra observation of fast variations of the 
hard X-ray emission in some spots in the remnant, possi- 
bly showing th e localized existence of even mG magnetic 
field strengths (lUchivama et al. . 2008f) (see, however also 
iButt et all <|2008h ). 



2 This is approximately also true if cross-field diffusion of 
the highest-energy nuclear particles modifies the shock almost 
everywhere, as discussed in the Introduction. 
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Fig. 1. Spatially i ntegrat ed, 7-ray spectral energy distribution of RX J1713. 7-3946 . The latest H.E.S.S. 7-ray data 
( Aharonian et all l2007al) ( blue color) are shown together with the EGRET upper limit for the RX J1713. 7-3946 
position (jAharonian et all l2006h (green color). The solid curve at energies above 10 7 eV corresponds to 7r°-decay 
7-ray emission, whereas the dashed and dash- dotted curves indicate the Inverse Comp ton (IC) and Nontherm al 
Bremsstrahlung (NB) emissions, respectively, from the theoretical model of BV06. See also Berezhko fc Volk ( 20071 ) . 



4. Comparison of the latest overall nonthermal 
spectrum with the BV06 spectrum 

We present in Fig. 2 along with new HESS data the latest 
hard X-ray spectrum, obtain ed with the Suzaku instru- 



ment ([Takahashi et al.l . 120081) . which is given in the form 



of an energy flux density bv luchivama et ah ( 20081 ). and 
compare it with the theore t ical sp ectrum of BV06 (see also 
Zirakashvili <k Aharonianl 12007 ). The Suzaku measure- 
ment does not cover the entire SNR, and therefore it needs 
to be suitably normalized by the require ment of optimum 



agreem ent with the ASCA spectrum, cf. lAharonian et al 



(2006), over the latter instrument's range between 0.5 and 
10 keV. The result is shown in Fig. 2. The good agreement 
basically stems from the fact that the amplified down- 
stream field of w 100 /iG, used to calculate the theoretical 
spectrum in BV06, already leads to electron synchrotron 
cooling above an electron momentum of ~ 500 m p c, and 
thus to a hardening of the synchrotron spectrum com- 
pared to a spectrum calculated without electron cooling 
(see Fig. 3 below). 

The same observed spectrum can also be compared 
with a theoretical spectrum (Fig. 3) in which a very low 
ion injection rate (77 = 10~ 5 ) and a rather low down- 
stream magnetic field strength of 20 fiG was assumed (see 
BV06). This corresponds to a dominantly leptonic 7-ray 



test particle spectrum without field amplification^]. The 
IC-scattered diffuse radiation field is the CMB plus inter- 
stellar far infra r ed and optical radiation fields as given in 
Berezhko et al.l (|2003bT ). This corresponds t o typical val- 



ues used for the Solar neighborhood ( e .g. iDrurv et al 
19941 iGaisser et all Il998t IPorter et all |2006|). We note 
that for nearby objects at distances d ~ 1 kpc the 
CMB contribution dominates in the IC emission (see also 
Porter et all l2006f) . 

The electron injection strength was fitted such that an 
optimum fit to the observations in the radio and X-ray 
ranges is achieved, cf. Fig. 3. We note that the form of 
the X-ray spectrum is only very poorly fitted in this lep- 
tonic scenario, especially when one uses the recent Suzaku 
measurements. Also the 7-ray spectrum has a maximum 
which is much too sharp in comparison with the observed 
H.E.S.S. spect rum. We note that, co mpared to earlier 
measurements ( Aharonian et alll2005l). the latest version 
of the H.E.S.S. 7-ray spectrum ([Aharonian et all 2007al ) 
deviates more clearly from the IC-type spectrum with a 
relatively sharp peak at e 7 ~ 1 TeV. Note also that our 



3 In fact, the strength of the downstream magnetic field 
might be even smaller by a factor of two or more, wherever 
the shock is not locally parallel to the external field. However, 
adopting such a small field would imply that even the gross 
amplitude of the maximum of the observed 7-ray energy flux 
could not be fitted to the observations. 
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Fig. 2. S patially integrated, over all nonthermal spectral energy distri bution of RX J17 1 3. 7-39 46 The ATCA radio 
data (cf. Aharonian et all |2006) [violet color], ASCA X-ray da ta (cf. Aharonian et all l2006l h Suzaku X-ray data 
' 2008! ) [red color], and 2006 H.E.S.S.-7 -rav data dAharonian et all l2007al) [blue color] are shown. The 



(|Uchivama et ah . ._ . _ i . 

EGRET upper limit for the RX J1713. 7-3946 position ( Aharonian et alT 2006h [green color] is shown as well. The solid 



curve at energies above 10 7 eV corresponds to 7r°-decay 7-ray emission, whereas the dashed and dash-dotted curves 
indicate the Inverse Compton (IC ) and Nonthermal Bremsstrahlung (NB) emissions, respectively, from the theoretical 
model of lBerezhko fc Volkl (|2006l ). 



spectrum of the nonthermal emission, that corresponds 
to the leptonic (or inefficient) scen ario, is almo s t iden - 
tical to the spectrum presented by IPorter et al. ( 2006I ). 
even though they approximate the electron spectrum an- 
alytically, whereas we calculate it numerically. Therefore 
it is also clear th at the quality of the fit achieved by 
IPorter et al. I (|2006h with a leptonic model will be consider- 
ably lower if one uses the Suzaku X-ray spectrum and the 
new HESS 7-ray spectrum instead of older, less accurate 
data. 

At 7-ray energies of 1 GeV the spectral energy flux 
density is a factor of about 30 below the value in the 
hadronic scenario. It might be possible to construct a 
more or less plausible form of the diffuse radiation field 
spectrum to obtain a reasonable fit in the TeV region. 
However, it remains very doubtful in our view whether 
this can avoid the sharp decline towards the GeV energy 
range indicated in Fig. 3. In any case such a construction 
cannot improve the unacceptable fit in the hard X-ray 
range. 

5. Conclusions 

We conclude that a theoretical acceleration model which 
takes into account magnetic field amplification and a con- 
sistent nuclear energetic particle production is consistent 
with the latest H.E.S.S. 7-ray and Suzaku hard X-ray ob- 



servations. It is also expected to be consistent with the ob- 
served good correlation between X-ray synchrotron emis- 
sion and VHE 7-ray emission. In the face of existing gas 
density variations in or near the SNR, a purely leptonic 
interpretation of the spatial correlation becomes poorer as 
this correlation becomes closer empirically. The attempt 
to explain the 7-ray emission by leptonic processes cannot 
be made consistent with the observed synchrotron spec- 
trum. The availability of higher-quality hard X-ray and 
7-ray measurements leads to a better consistency with 
the kinetic nonlinear theory prediction, whereas the con- 
sistency of the leptonic model becomes much poorer. 

Whether a leptonic model can lead to an acceptable fit 
of the 7-ray spectrum through a more detailed evaluation 
of the local diffuse radiation field is an open question. In 
our view, even such a modification would encounter enor- 
mous difficulties in attempting to fit the 7-ray spectrum 
over the additional two orders in 7-ray energy down from 
the VHE range to the GeV range. The leptonic scenario 
seems also to be inconsistent with the filamentary X-ray 
morphology which suggests substantial field amplification 
at least over part of the remnant. 
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